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A review of the calculation of the two loop bosonic corrections to Ar is presented. Factorization and matching 
onto the Fermi model are discussed. An approximate formula, describing the quantity over the interesting range 
of Higgs boson mass values from 100 GeV to 1 TeV is given. 


The muon decay lifetime (r^) has been used 
for long as an input parameter for high precision 
predictions of the Standard Model (SM). It allows 
for an indirect determination of the mass of the 
W boson (Mw), which suffers currently from a 
large experimental error of 33 MeV [1], one or¬ 
der of magnitude worse than that of the Z boson 
mass {Mz)- A reduction of this error by LHC 
to 15 MeV [2] and by a future linear collider to 
6 MeV [3] would provide a stringent test of the 
SM by confronting the theoretical prediction with 
the experimental value. 

The extraction of M\y with an accuracy match¬ 
ing that of next experiments, i.e. at the level of 
a few MeV necessitates radiative corrections be¬ 
yond one loop order. Large two loop contribu¬ 
tions from fermionic loops have been calculated 
in [4]. The current prediction is affected by two 
types of uncertainties. First, apart from the still 
unknown Higgs boson mass, two input parame¬ 
ters introduce large errors. The current knowl¬ 
edge of the top quark mass results in an error 
of about 30 MeV [5], which should be reduced 
by LHC to 10 MeV and by a linear collider even 
down to 1.2 MeV. The inaccuracy of the knowl¬ 
edge of the running of the fine structure constant 
up to the Mz scale, Aa{Mz), introduces a fur¬ 
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ther 6.5 MeV error. Second, several higher order 
corrections are unknown. In fact the last cor- 
rection of order O(a^), generated by purely elec¬ 
troweak bosonic loops, has been calculated only 
recently in [6] and confirmed by an independent 
group in [7]. The details of both calculations are 
presented in [8]. 

The two loop bosonic corrections to the muon 
lifetime have been previously estimated by two 
different methods. First, the leading term in the 
large Higgs boson mass expansion has been ob¬ 
tained in [9] and [10]. Although the results were 
different, the size of the contributions was neg¬ 
ligible. Apart from the disagreement, the large 
mass expansion is not justified for small Higgs 
boson masses, and the experimental data seems 
to favour a low mass range. Second, a resumma¬ 
tion formula has been used [5], which with respect 
to bosonic corrections consists of assuming a ge¬ 
ometric progression of the contributions with the 
order of perturbation theory. This means that 
at two loop order, the bosonic corrections should 
be equal to the square of the one loop result. 
This, however, is contradicted by the fact, that 
the leading behaviour with the Higgs boson mass 
would then be of logarithmic type log^(M^)), 
whereas it is known that the result should behave 
as the square of this mass [11]. 

Since both estimates turn out to be unreliable 
the full calculation was necessary. Here we shall 
describe the crucial ingredients, starting from the 
proper definition of A?’, passing through the var- 
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ious stages of the renormalization procedure, and 
ending with the methods used in numerical eval¬ 
uation. 

The muon decay is described by an effective 
field theory, the Fermi model, which consists of 
QED, QCD of five light quarks and the four 
fermion point interaction 


Gf_ 

-Cj’=-^e7“(l - 75)/x(g) r'^7a(l - 75)i^e- (1) 


The constant Gp should be predicted by the 
Standard Model. The difference between the tree 
level prediction and the full result is factorized in 
the quantity Ar 


V2 


7ra 
2 m2 


(1 -f Ar). 


( 2 ) 


In this sense Ar can be considered as the match¬ 
ing or Wilson coefficient in the Fermi model. The 
matching equation assumes the form 


^SM ^ g 
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( 3 ) 


where p and m denote the external momenta and 
the light masses respectively, and S is the renor¬ 
malized scattering matrix in both models. The 
easiest way to obtain Ar is to simply put p and 
m equal to zero on both sides of the equation. 
This procedure will of course generate spurious 
infrared divergences, which however will cancel 
from the equation as has been shown in [12]. The 
main advantage comes from the fact, that on the 
right side, all loop diagrams will be scaleless (as 
they contain only light particles), and will there¬ 
fore vanish, leaving only the tree level diagram 
trivially proportional to Gf- Additionally, the 
SM diagrams on the left hand side will reduce to 
vacuum bubbles. 

A major point in the above procedure is con¬ 
nected with the treatement of spinor chains. Di¬ 
mensional regularization for spurious infrared di¬ 
vergencies forbids a direct use of the Chisholm 
identity 


= gpulp + 9uplp - gpplu - iepupcrl'^iy) 

for all diagrams, which are divergent after renor¬ 
malization. A priori we would have to introduce 



Figure 1. Example of a Fierz rearrangement 
about the photon line. 


a whole basis of tensor products of possible spinor 
chains. It turns out, however, that the projection 
onto the Fermi operator can be defined unam¬ 
biguously through Fierz symmetry. The problem¬ 
atic diagrams contain always a photon line con¬ 
necting the two chains. A Fierz rearrangement 
about this line or the W boson line, as depicted 
in Fig. I, leads to an object which has the re¬ 
quired vertex structure 

I^Pl^IpPl- (5) 

In practice, this Fierz rearrangement can be re¬ 
alised with the help of a projector. In the present 
case, the transformation 

® ^ ~ 2D{D 

gives directly the coefficient of the operator Eq. 5 
in D dimensions. 

The matching procedure requires proper renor¬ 
malization of the model. Although in the end we 
are interested in the on-shell parametrization for 
the masses and the electric charge, an intermedi¬ 
ate renormalization can be chosen at will. In fact 
the factorization theorem [12] has been proven 
for mass independent schemes and the on-shell 
scheme does not have this property. However, as 
long as the parameters are defined at the heavy 
scale the mass dependence does not contradict the 
theorem. For this reason, we may safely renor¬ 
malize the boson masses in the on-shell scheme. 
The situation is somewhat more problematic with 
the electric charge, which is usually dehned by the 
Thompson scattering processes for on shell elec¬ 
trons and photon. If we are interested in bosonic 
corrections only, the definition can be kept, since 
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an inspection of the diagrams leads to the con¬ 
clusion that no spurious infrared divergencies can 
be introduced. For light fermionic contributions 
the situation is different with the usual solution 
consiting of a shift of the definition to the Mz 
scale [13]. 

In this work, two renormalisation schemes have 
been used. First, the complete calculation has 
been performed in the on-shell scheme. Sec¬ 
ond, the model has been renormalized in the MS 
scheme, and the result translated back to the on- 
shell scheme by means of relations between the 
on-shell and MS masses, which for the W and 
Z masses were required up to two loop order. 
Obviously the translation should be applied to a 
scheme independent quantity. In this case, this is 
the Fermi constant Gp, and the correct relation 
is 


2M2 


(1 -I- Ar) 


7ra 


OS 


2M2 


(1 + Ar) 


■( 7 ) 
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It is not trivial that the formulae will coincide, 
since the on-shell counterterms on the left hand 
side contain terms of order e (e = 1/2(4 — ZI)), 
which are not used in the translation between the 
schemes on the right hand side. 

An interesting part of the definition of the 
model is the treatement of tadpole diagrams. It 
is known that gauge invariance of mass counter¬ 
terms requires inclusion of tadpoles [14,15] (at the 
two loop level this has been explicitely shown in 
[16]). In this case, however, one cannot use one- 
particle-irreducible (IPI) Green functions. In or¬ 
der to have gauge invariant counter-terms and 
IPI Green functions only, a special procedure 
was designed. An additional renormalisation con¬ 
stant for the bare vacuum expectation value uq, 
denoted has been introduced and explicitely 
split from the bare masses 


vq —> vqZI/'^, ( 8 ) 

^W,Z ( 9 ) 

The term linear in the Higgs field H in the la- 
grangian 

T°H° = - l)H°, (10) 

eo 


is then used to determine through the require¬ 
ment that tadpoles are canceled. It can be proved 


[8,14] that the bare masses are gauge invariant in 
this case. 

It should be stressed that this procedure can 
be advantageous at higher orders, since the num¬ 
ber of diagrams is strongly reduced by using IPI 
Green functions. For the muon decay, there are 
twice as many diagrams if tadpoles are intro¬ 
duced. If a treatement per diagram is required 
{e.g. expansion), the calculation time is strongly 
correlated to the number of diagrams and sub¬ 
stantial differences in execution can be observed. 

Having defined the method of factorization and 
renormalization, we turn to the calculation of the 
respective diagrams. As noticed already above, 
all of the bare diagrams are of the vacuum type, 
and for these a reduction procedure supplied with 
analytical formulae is known [17]. The situation 
is slightly more complicated with the on-shell two 
point functions that are needed for the on-shell 
definition of masses. The strategy adopted in this 
work consists of two algorithms. The first one [18] 
reduces the tensor integrals to scalar ones and 
uses topological symmetries to lead to as small a 
set of basic integrals as possible. The integrals 
obtained in this way are now subject to numeri¬ 
cal evaluation. We used the efficient one dimen¬ 
sional representation given in [19] and extended 
its implementation (S2LSE) to work with soft¬ 
ware emulated quadruple precision. 

Several tests have been applied to the numer¬ 
ics. First of all, the Z boson propagator has been 
completely tested by means of the low momentum 
expansion, which is valid in this case, as all dia¬ 
grams are below threshold. A 6 digit coincidence 
has been found at tenth order. The W propa¬ 
gator required the use of large mass expansions, 
since several diagrams are at threshold, but the 
same precision as previously was also obtainable. 
Finally both propagators were tested against the 
result for the MS to on-shell mass relations of 
the two gauge bosons given in [16]. This was also 
considered as a test of the gauge invariance of the 
mass counterterms. 

Concerning gauge invariance, all of the calcu¬ 
lations apart from the two two loop mass coun¬ 
terterms, have been performed in the general 
gauge with three parameters and the cancellation 
of the dependence on these parameters has been 
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Figure 2. The full result for the bosonic correc¬ 
tions to Ar (solid line) against the large mass ex¬ 
pansion with leading term and the square of the 
one loop result. 


explicitely verified. 

The full result is given in Fig. 2 as function of 
the Higgs boson mass, for the following values of 
the input parameters 

0-1 = 137.03599976, (11) 

Mw = 80.423 GeV, 

Mz = 91.1876 GeV. 


On the same plot, the leading term of the large 
Higgs boson mass expansion [6,9] and the square 
of the one loop result are shown to contrast the 
previous estimates with the full result. 

Our result has also been reexpanded in the 
large Higgs boson mass for comparison with [7] 
using the MS value and the expansion of the 
translation formulae given in [16]. The five term 
expansion is also shown in Fig. 2. 

It turns out that the variation of the result with 
the W mass within the experimental error bars of 
±33 MeV [1] is negligible. The correction can be 
described to 5% accuracy outside of the double 
threshold region for Higgs boson masses ranging 
from 100 GeV to 1 TeV, with the following for- 


Figure 3. W boson mass shift generated by two 
loop bosonic corrections. 


mula 

d ± 5A ± ciT ± C 2 T , 


where 


A = 


Mh 

100 GeV 


2 

- 1 , 


L = log 


/ Mb y 
V100 GeV J ’ 


and the coefficients are 

a = 6.385 X 10-^ 

6 = -1.110x10“®, 
Cl = -2.870 X 10“®, 
C2 = 6.445 X 10“®. 


( 12 ) 


(13) 

(14) 


(15) 


The effect on the prediction of the W boson 
mass from the muon decay lifetime is small and 
can be obtained from an expansion. In fact 6% 
accuracy is guaranteed by estimating the result¬ 
ing mass shift from the formula 

AMw = -1.491 X 10'‘Ary^ [MeV]. (16) 

The respective plot is shown in Fig. 3. 
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In conclusion, the calculation of the two loop 
bosonic corrections to the muon decay lifetime 
as parametrized by Ar has been reviewed. The 
crucial steps of the factorization and matching 
have been described together. A simple formula 
describing the correction for a large range of Higgs 
boson masses has been given. 
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